Likelihood Approach to the First Dark Matter Results from XENON100 
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Many experiments that aim at the direct detection of Dark Matter are able to distinguish a 
dominant background from the expected feeble signals, based on some measured discrimination pa- 
rameter. We develop a statistical model for such experiments using the Profile Likelihood ratio as a 
test statistic in a frequentist approach. We take data from calibrations as control measurements for 
signal and background, and the method allows the inclusion of data from Monte Carlo simulations. 
Systematic detector uncertainties, such as uncertainties in the energy scale, as well as astrophysical 
uncertainties, are included in the model. The statistical model can be used to either set an exclusion 
limit or to quantify a discovery claim, and the results are derived with a proper treatment of statis- 
tical and systematic uncertainties. We apply the model to the first data release of the XENON100 
experiment, which allows to extract additional information from the data, and place stronger limits 
on the spin-independent elastic WIMP-nucleon scattering cross-section. In particular, we derive a 
single limit, including all relevant systematic uncertainties, with a minimum of 2.4 x 10 -44 cm 2 for 
WIMPs with a mass of 50GeV/c 2 . 
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I. INTRODUCTION 



It is now well established that the dominant mass frac- 
tion of our Universe consists of some yet-unknown form 
of dark matter [1]. Well- motivated models predict Dark 
Matter in the form of Weakly Interacting Massive Parti- 
cles (WIMPs) [2] for which searches can be conducted in 
direct scattering experiments located on Earth. In such 
experiments, WIMPs are expected to induce nuclear re- 
coils with a roughly exponentially falling spectrum that 
extends to at most a few tens of keV in energy [3]. In 
contrast, the dominant background is usually composed 
of electronic recoils from f3- or 7-radiation. Most direct 
Dark Matter search experiments are able to discriminate 
nuclear recoils from electronic recoils based on some dis- 
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crimination parameter, such as the ratio of ionization to 
phonon signal [4, 5], scintillation to phonon signal [6], 
ionization to scintillation signal [7, 8], or the pulse shape 
of the scintillation light [9]. In a plot of this discrimi- 
nation parameter as function of energy, the signal and 
background events thus separate. Typically, a signal ac- 
ceptance area is defined a priori in this parameter space. 
Only events that fall within this area are considered as 
WIMP candidate events for the interpretation of the re- 
sult in terms of a Dark Matter detection, or for the cal- 
culation of limits on the WIMP-nucleon scattering cross- 
section. 

Such a hard cut has some obvious disadvantages. Most 
importantly, the particular location of the signal can- 
didate events within the signal acceptance area is not 
taken into account. In general, the expected signal and 
background will follow a different spectrum, which is uti- 
lized by the maximum-gap and optimum-interval meth- 
ods [10], a frequentist approach to calculate limits. The 
generalization to two dimensions, the maximum-patch 
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method [11], extends this feature to the case where back- 
ground leaks into the signal acceptance region from high 
or low discrimination parameters. An intrinsic advan- 
tage of these methods is that the resulting limit is ro- 
bust against small changes of the signal acceptance re- 
gion, and they are widely used in discriminating direct 
Dark Matter searches today. However, these methods 
also have some severe drawbacks: they are designed to 
always result in an upper limit, and a natural transition 
to a claim of detection, as for example inherent to the 
Fcldman-Cousins-approach [12], is not possible. In addi- 
tion, systematic uncertainties are not taken into account. 
Thus, the resulting limit can be quoted at a given statis- 
tical confidence level (CL), typically taken to be 90%, but 
systematic uncertainties need to be treated separately. 

In this paper, we present an approach that overcomes 
both of these problems of the conventional treatment of 
data from direct Dark Matter detection experiments. As 
a practical example, we use data from the first Dark 
Matter results of the XENON100 experiment [8]. This 
data was previously analyzed in the classical way, by pre- 
defining a signal acceptance region based on calibration 
data. In particular, this region was defined using the 
prompt scintillation signal as a measurement of recoil en- 
ergy, in the range between 4-20 photoelectrons (PE). In 
energy, this corresponds to about 9 — 33 keV nr (keV nu- 
clear recoil) , depending on the energy calibration through 
the assumed relative scintillation efficiency £ e s [13, 14]. 
The logarithm of the ratio of proportional scintillation, 
S2, to prompt scintillation, SI, was taken as a discrim- 
ination parameter [15], and the signal acceptance region 
was defined below the median of the nuclear recoil band, 
as measured with a AmBe neutron calibration. The 
electronic recoil discrimination in this signal acceptance 
region was determined to be larger than 99% based on 
calibration data from Compton-scattered 60 Co gammas. 
Using 11.17 live days of data, no events were observed 
in the signal acceptance region. Thus, upper limits on 
the WIMP-nucleon scattering cross-section were derived, 
based on zero observed events. Two separate limit curves 
were calculated for two different parameterizations of C c g 
in the low-energy region where measurements are lacking 
or are uncertain [16]. 

Here, we re-analyze this data with the Profile Likeli- 
hood Ratio statistical approach [17], that is inherent to 
the MINOS package commonly used with the MINUIT 
software [18]. This method relies on the assumption that 
one can properly model the background via calibration 
measurements or Monte Carlo simulations, as opposed 
to the Yellin methods [10], where no information on the 
background is assumed. This increases the sensitivity 
of the method at the price of being less robust against 
systematic uncertainties on the background. However, 
by incorporating systematic uncertainties into the model, 
one can adjust the trade-off between sensitivity and ro- 
bustness to the level desired by the experimentalist. In 
any case, understanding of the background is of course 
essential to any method that is used to quantify the sig- 



nificance of a discovery claim. 

The particular Profile Likelihood model used here is 
presented in Section II. Uncertainties in C e g, that can 
serve as an example for systematic uncertainties related 
to the detector, are treated in Section III, together with 
systematic uncertainties from the escape velocity t> CS c as 
an example for astrophysical uncertainties. Calibration 
data for both the electronic recoil background as well as 
the nuclear recoil signal are taken as constraining control 
measurements. These constrain the full likelihood model, 
which is constructed in Section IV. The sensitivity of the 
experiment is calculated, based solely on this model, in 
Section V. Taking the actual measurement into account, 
the resulting exclusion curve is derived also in Section V. 
Discussion and conclusions follow in Section VI. 



II. THE STATISTICAL MODEL 

Events recorded by the XENON100 detector can be 
characterized by their prompt and proportional scintil- 
lation signals, and are therefore treated in this analysis 
as a set of points in the (SI, S2) plane. The distribution 
of signal events in SI can be predicted from theoretical 
models taking into account the detector response, as will 
be described below. The joint distribution in 51 and S2, 
for signal and background events, is estimated from the 
calibration measurements. The current measurements of 
£ e ff and v esc are similarly treated as another set of con- 
trol measurements. The statistical model describing the 
outcome of the measurements thus includes a number of 
unknown quantities, i.e. nuisance parameters. Those are 
the expected number of background events Nj, given the 
exposure and the (SI, S2) parameter range in question, 
a set of probabilities e s = {e^} and e b = {e^} describing 
the distribution of signal and background events in the 
(Si, S2) plane, the relative scintillation efficiency C e s, as 
well as the escape velocity v csc - The full likelihood func- 
tion , for a given WIMP mass m x and cross-section a, 
is written as a product of five terms: 

Jz? = J>?i(cr,N b ,e s ,e b ,£ cS ,v csc ;m x ) (1) 
xJf 2 (e s ) x Jzf 3 (e b ) 
xif 4 (£ off ) x _Sf 5 (u esc ). 

The first term Jzfi describes the main measurement of the 
XENON100 detector, while the following terms describe 
the subsidiary measurements that are used to constrain 
the nuisance parameters in the main likelihood term Jz?i . 
The precise definition of these terms will be given in the 
following sections. Additional uncertainties can easily be 
incorporated by adding additional Likelihood terms, but 
are not relevant for the data set analyzed here. 

The signal cross-section a is the one parameter of in- 
terest. All other parameters are nuisance parameters 
which are profiled out with a profile likelihood ratio, as 
explained below (Equation 3). We follow the procedure 



3 



of a hypothesis test based on the profile likelihood ra- 
tio [19, 20]. This technique can be used both to exclude 
a WIMP with a specific mass and cross-section, or to 
establish the significance of a discovery. 



A. Exclusion 

A test statistic q a reduces the observed data to only 
one value and is constructed in order to test the signal 
hypothesis H a . It is given by 



-2 In A(er) a < a 
a > a 



(2) 



where a is the maximum likelihood estimator (MLE) of 
tr, i.e. the value of a that maximizes the likelihood Equa- 
tion 1. A(cr) is the Profile Likelihood ratio and is given 
by 



A(<t) 



max Jif (a;£ cS ,v csc , N b ,e s ,e b ) 

(7 fixed 

max (tr, C c « , w osc , N b ,e s ,e b ) 
[a, C Q g , Wcsc, N b , e s , e b 



Ve SC ,N b , e s ,i b 



(3) 



The double-hat parameters in the numerator are the con- 
ditional MLEs of the nuisance parameters when the sig- 
nal cross-section is fixed to a given value a. The 'single- 
hat' parameters in the denominator are the MLEs of all 
parameters allowing also er to vary. By construction, 
< A(cr) < 1, hence q a > 0. q a equals zero when the 
best-fit value of the cross-section (a) equals the hypoth- 
esized value (er), which corresponds to the most signal- 
like outcome. Larger values of the test statistic q a in- 
dicate that the data are less compatible with the signal 
hypothesis H a . Since we are concerned with calculating a 
one-sided upper bound, we only consider outcomes with 
a < a as an evidence against the signal hypothesis and 
set q a to zero otherwise. 

Let f{q a \H a ) be the probability distribution function 
of the test statistic q a under the signal hypothesis H a , 
and let g° bs be the value of the test statistic obtained 
with the observed data. The signal p- value p s , is the 
probability that the outcome of a hypothetical, random 
XENON100 experiment results in a test statistic larger 
(less signal-like) than the observed one, when the signal 
hypothesis H a is true. Therefore, p s given by 



Ps = I f(q*\H,j)dq 



(4) 



The signal hypothesis H a is rejected at 90% CL if p s < 
10%. 

Downward fluctuations of the background might lead 
to exclusions of very small cross-sections to which the 



experiment is not sensitive. To protect against such an 
effect, p s is modified [21, 22] to 



P, 



Ps 



I-Pb 



where 



1 ~Pb = I f(qa\H )dq a 

'« b! 



(5) 



(6) 



is the probability of the test statistic q a to be larger than 
the observed one under the background-only hypothesis 
H a= Q = Hq. Other protection procedures would also 
be conceivable, but this particular extension has a con- 
servative over-coverage nature. It has been verified by 
Monte Carlo simulations that the coverage of our claimed 
90% CL upper bound on the cross-section is in the range 
(92-95)% and thus larger than 90%, in accordance with 
the conservative nature of the method. 

The upper limit a up (m x ) on the cross-section a for a 
given WIMP mass m x is found by solving 



P> 



r up 



{m x )) = 10%. 



(7) 



Wilks' theorem [23] states that q a follows a chi-square 
distribution in the limit of a large number of observa- 
tions, which in this context not only includes the ob- 
served WIMP candidate events, but also all of the control 
measurements. It has been verified by Monte Carlo sim- 
ulations that this asymptotic behavior indeed describes 
the distribution of the test statistic to a good approxi- 
mation. We therefore use the chi-square approximation 
in order to estimate the signal p- value p s . p b is estimated 
from Monte Carlo simulations of background only events. 



B. Discovery 

We can also use the above statistical method in a nat- 
ural way to quantify the significance of a possible addi- 
tional event population or signal discovery. To this end, 
we test the background-only hypothesis (a = 0), and try 
to reject it. Similar to Equation 2, the discovery test 
statistic go is defined as 



9o 



-21nA(0) <7>0 
a< 0. 



(8) 



Setting (7 = in Equation 3 yields 

Jzf ( a = 0, £ e ff, "esc, Nb, e s ,e b 

A(0) 



(a, C c s, Vcsc, N b , l s ,h 
The p- value is defined as 



(9) 



Po = I f(qo\H )dq , 

/ g obs 



(10) 
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where f(qo\Ho) is the probability density function of qq 
under the background-only hypothesis Hq. A p- value of 
Po = 0.00135 corresponds to a 3a signal evidence. Fol- 
lowing Wilks' theorem, go asymptotically distributes as 
a chi-square distribution under the background-only hy- 
pothesis. Under Wilks' theorem, the discovery signifi- 
cance is given by Za with Z = \f q obs . As a general 
remark, the Profile Likelihood can result in the discovery 
of a process that looks like the signal hypothesis. How- 
ever, the interpretation of the physical origin of such a 
signal is of course left to scientific discourse. 



III. EXPECTED SIGNAL WITH £ cff AND v esc AS 
CONTROL MEASUREMENTS 



with i bs = representing the £ e ff-median in Figure 1. 
t = — 1 corresponds to the lower boundary of the 68% 
confidence region of C e g (dark shaded band in Figure 1), 
etc.. 

The number N s of expected signal events is given by 
the integral over the predicted WIMP energy spectrum 
dR/dE nr [34] as seen in the detector. dR/dE DI is calcu- 
lated for a given WIMP mass m x , cross-section a, Galac- 
tic escape velocity i> esc as well as a given target mass and 
exposure time. 

For the uncertainty on the escape velocity v csc , we use 
the asymmetric probability distribution /„ as shown in 
Figure 7 of [35], yielding 498km/s < i; GSC < 608km/s at 
90% CL, with the median being v ^ s — 544km/s. The 
likelihood term is defined accordingly to be 



The relative scintillation efficiency £ c g is the largest 
systematic uncertainty for the data analyzed here. It 
relates the expected number of SI photoelectrons to the 
recoil energy E nl -. The recoil energy dependence of £ e ff 
together with its uncertainty are taken from a Gaussian 
fit of all available direct measurements [24-30] as shown 
in figure 1 . To extrapolate below 5 keV nr , where there are 
no available measurements of £ g yet, the 68% confidence 
interval for C c g is constructed such that it follows the 
two different extrapolations used in [8] to facilitate the 
comparison of the resulting limit. We emphasize that 
this extrapolation follows both the new data published 
in [31], the theoretical expectations spelled out in [32], 
as well as the trend expected from simulations [33] . 
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■&b{Vesc) = fv{v bs\ V e 



(12) 
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FIG. 1: All direct measurements of C e s [24-30] together with 
a Gaussian fit (solid blue line) and the uncertainty band indi- 
cating the ±1 and ±2<r regions (shaded dark and light blue, 
respectively). 

We parametrize the uncertainty of £ c g with a single 
nuisance parameter t, which is defined to be normally 
distributed with zero mean and unit variance. Thus, the 
£ e ff-term in Equation 1 corresponds simply to 



Let v(E nr ) be the expected number of photoelec- 
trons (PE) for a given recoil energy E nT , given by 

S 

u(E nr ) = E m x £ cS (E nr ) x x L y . (13) 

*^ee 

Sec = 0.58 and S nl — 0.95 [15] are the scintillation 
quenching factors due to the electric field for electronic 
and nuclear recoils, respectively. The normalization 
light yield L y for this data sample was measured to be 
L y (122keV cc ) = (2.20 ± 0.09) PE/keVec [8]. Errors on 
these parameters are small and have been neglected for 
the current analysis. 

The signal rate in number of photoelectrons n is then 
given by 



dR 
dn 



dE — (E;m x ,v csc ) x Poiss [n\v{E)) . (14) 



Taking into account also the finite average single- 
photoelectron resolution ctpmt = 0.5 PE of the 
XENON100 photomultipliers, the resulting Sl-spectrum 
is given by 



dR 
d~Sl 



> Gauss(5*l|n, \Aio"pmt) 



n=l 

dR 

x 1 x Ccuts 

dn 



(SI) 



(15) 



where Ccuts(Sl) is the acceptance of the applied cuts [8]. 

The total number of expected signal events 
N s (t,v csc ,a) is calculated from the integral 



N s (t,v e 



Sly. 

Silo, 



dR 
dSl 



(m x ) dSl 



(16) 



jSf 4 (£eff(u) = <op(-(t - iobs) 2 /2) 



(11) 



where Slower = 4 PE and iSl uppcr = 20 PE is the energy 
interval considered in the analysis. The effect of varying 
t (Equation 11) for different WIMP masses is shown in 
Figure 2. It is calculated for an exposure of 11.17 days 
with 40 kg of xenon and a cross-section of a = 10 -39 cm 2 . 
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FIG. 2: Expected number of signal events N a as a function of 
the £ e ff parametrization t, for various WIMP masses m x as 
indicated in the figure, and for a cross-section a = 10~ 39 cm 2 . 
The lower the WIMP mass m x , the stronger the systematic 
uncertainty from C c s. 



Finally, the predicted normalized WIMP spectrum 
/ s (51) is given by 



/ s (51;m x ) 



dR 
dSl 



(m x ) 



si a 



Sh, 



dR 
dSl 



(17) 



(m x ) d51 



Variation of t modifies the expected number of total sig- 
nal events at a given cross-section, but the shape of the 
51 distribution is found to be hardly affected, see Fig- 
ure 3. We therefore assume that f s is independent of t, 
which greatly simplifies the analysis. 
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FIG. 3: The normalized WIMP recoil spectrum f s (Sl) for 
m x — 5GeV/c 2 and three possible parameterizations of £ e ff 
through the parameter t (Equation 11). Even for low WIMP 
masses, the spectral shape hardly changes. 



IV. LIKELIHOOD CONSTRUCTION 

The total number of calibration events used for this 
analysis in the 4 PE-20 PE 51-region was M s = 18907 
nuclear recoils from 241 AmBe and Mb — 5053 electronic 
recoils from Compton-scattered 60 Co gammas. Those 
were divided into j = 1, . . . , (K = 23) bands in the 
(51,52) plane, containing approximately equal fractions 
(about 4%) of nuclear recoils in each band, see Figure 4. 
e J s and t\ are the corresponding probabilities for a signal 
or background event to fall in a given band j. These 
bands were constructed such that the fraction of nuclear 
recoil events in each band is independent of 51, so that 
the signal 51 spectrum f s is the same in all bands. The 
binning resolution of the (51, 52)-plane is limited by the 
amount of available neutron calibration data. It was veri- 
fied that the resulting limit is not sensitive to the number 
of bands. 
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FIG. 4: Scatter plot of electronic recoils from a 60 Co source 
(blue circles) and nuclear recoils from an 241 AmBe source (red 
dots). Superimposed are the border lines of the bands j (thin 
lines) as well as the median of the electronic (upper thick line) 
and nuclear (lower thick line) recoil bands in the relevant 5*1 
parameter range (dashed vertical lines). 

We take the calibration measurements as control mea- 
surements of m? s and m? h events in band j, in order to 
constrain e 3 s and t\. 



and 



are assumed to be Pois- 



son distributed with expectations e 3 s M s and e J h Mb. The 
corresponding likelihood terms in Equation 1 are thus 

K 

= \[\><»»:n>:<l.\Ll (18) 

3 

K 

Mi) = n Poiss KK M &)- ( 19 ) 



Finally, given a set of n 3 data points (51j,52j) in a 
band j, the likelihood that they emerge from a given 
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WIMP spectrum is given by 

K 



JSfi - Y[Poiss(^\eiN s + eiN b ) 

3 



1=1 



eiN s + e 3 h N b 



where f s is given by Equation 17. In accordance with 
both observation and Monte Carlo simulations, the elec- 
tronic recoil background spectrum ft is assumed to be 
flat in energy [36] . 

Collecting all the likelihood terms together, the full 
likelihood function reads 



K 



5? = []Poiss(^|e^ s +e J b iV b ) 



n 



eiN s f s (Sl) + elN b f b (Sl) 



elN s + elN b 



xPoiss(m{\e J b M b ) x Poiss(m J s \e J s M s ) 



xe 



-(t-t obB ) 2 /2 



X fv { v obs l^esc) ■ 



(21) 



Differentiating logJ>f with respect to the expected num- 
ber of signal events N s and with respect to e J e , e J b , we find 
the following relations between the MLEs: 



K 

N = N S + N b = J2 nJ 

3=1 



(22) 



(N s + M s y4 + {N b + M b )e{ = n? + mi + ml (23) 

These relations prove useful in the minimization of the 
likelihood. 



V. WIMP EXCLUSION 

To assess the statistical power of the experiment, the 
exclusion sensitivity can be calculated even before ana- 
lyzing the actual data. The sensitivity of the experiment 
is the limit that can be expected for given exposure and 
experimental conditions. The test statistic q a is obtained 
by plugging the likelihood (Equation 21) into equations 2 
and 3. To verify the validity of Wilks' theorem, we gen- 
erate a signal for a given WIMP mass and add it to a 
background simulation. The resulting test statistic dis- 
tribution f(q a \Ha) for a 5GeV/c 2 WIMP is shown in 
Figure 5 (thick red histogram). One can clearly see that 
the distribution is well approximated by a chi-square dis- 
tribution. This allows to estimate p s analytically. 

To estimate the exclusion sensitivity, background-only 
experiments are simulated, based on a Poisson distribu- 
tion with the expectation of 22 events as observed in [8] . 
The test statistic distribution f{q a \Ho) is shown as the 
thin (blue) curve in Figure 5. The sensitivity is defined 




FIG. 5: The probability density functions f{qa\H a ) (thick red 
histogram) and f(q a \Ho) (thin blue histogram) for a WIMP 
with m x = 5GeV/c 2 and a = a up . One can clearly see the 
distribution is well approximated by a chi-square distribution 
(dashed red line). The median value of f(q a \Ho) is indicated 
together with the test statistic q& bB observed in data. 



as the median medf^] of f(q a \Ha), also indicated in 
the figure. One can see that a large fraction of exper- 
iments under the background-only hypothesis (thin blue 
histogram) result in a test statistic q a that is similar to 
or even more signal- like (the area to the left of q% bs ) than 
that of experiments with a signal (thick red histogram) . 
Figure 6 shows the sensitivity as its la (dark shaded) and 
2cr (light shaded) bands. For completeness of the statis- 
tical analysis we note that the p- value of the background- 
only hypothesis is 50% as there are less events than ex- 
pected, a < and qa = 0. 

To derive the upper bound on the WIMP cross-section 
as function of WIMP mass, we solve Equation 7 for the 
cross-section a = a up (m x ) that satisfies p' s — 10% (see 
Section II). The test statistic g° bs actually observed in 
data is shown in Figure 5, and the actual limit is shown 
in Figure 6. This limit has a minimum of cr up — 2A x 
10- 44 cm 2 at m x - 50GeV/c 2 . 

Two major advantages of the method presented here 
manifest themselves when comparing this new limit to 
the conventionally calculated one previously published 
in [8] . In the high WIMP mass region, the limit obtained 
here is stronger by about 20%-30%. This simply reflects 
the fact that in [8], only half of the available parame- 
ter space (50% nuclear recoil acceptance) was defined a 
priori as the signal acceptance region. In contrast, the 
current analysis considers the full available discrimina- 
tion parameter space, also taking the background in each 
band j into account. Since there are no events observed 
in the data even in the la region above the nuclear recoil 
median (see Figure 3 in [8]), the limit improves accord- 
ingly. 

The second major difference between these two analy- 
ses appears in the low WIMP mass region, which is most 
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FIG. 6: Parameter space of spin-independent elastic WIMP- 
nucleon cross-section a as function of WIMP mass m x . The 
sensitivity for the data set analyzed here is shown as light 
and dark (blue) shaded areas at la and 2a CL, respectively. 
The actual limit at 90% CL, taking into account all relevant 
systematic uncertainties as derived with the Profile Likelihood 
method, is shown as the thick (black) line. Two limits from 
the same data set, derived for two assumptions of the behavior 
of £ e ff , are shown as dotted lines [8]. A limit from CDMS [4] is 
shown as thin (orange) line, re-calculated assuming an escape 
velocity of 544km/s and Vo = 220km/s. Expectations from 
a theoretical model [37] are shown, as well as the areas (at 
90% CL) favored by CoGeNT (green) [38] and DAMA (red, 
without channeling) [39]. 



lation of the limit and results in a single limit at 90% CL, 
including all considered uncertainties. 



VI. CONCLUSIONS 

We have introduced and applied the frequentist statis- 
tical method based on the profile likelihood test statistic 
to re-analyze the first data release from the XENON100 
direct Dark Matter search experiment. This method 
avoids the need to a priori define a signal acceptance 
region, but instead takes all measured data into account. 
The background was estimated using calibration data 
as control measurements. Uncertainties in the relative 
scintillation efficiency C Q s and the Galactic escape ve- 
locity v osc were taken into account in the construction 
of the likelihood model. Using the profile likelihood test 
statistic allows to calculate the sensitivity of the exper- 
iment with its uncertainty bands, and to set a well- 
defined single limit, taking all systematic uncertainties 
into account. Applying the method to the previously 
published XENON100 data results in an improvement of 
the limit over a wide WIMP mass range, with a min- 
imum cr up < 2.4 x 10~ 44 cm 2 for WIMPs with mass 
rn x = 50GeV/c 2 . In addition, this method can easily 
be applied for the discovery of a WIMP signal. 



sensitive to uncertainties in £ e ff. The limit obtained here 
is more stringent than the limit shown in [8] , constraining 
the areas favored by CoGeNT [38] and DAMA [39, 40]. 
This comes from our definition of the uncertainty bands 
of £ e g as shown in Figure 1, where we set the ler in- 
terval to equal the two extrapolations used in the previ- 
ous analysis. We emphasize that here we do not assume 
any specific behavior of £ c g, but rather consider it as 
an unknown quantity with the appropriate measurement 
uncertainty, as described in Section III. Such systematic 
uncertainty is incorporated in a natural way in the calcu- 
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